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Abstract 

We solve in this paper the problem of the antiferromagnetic transition for the Q-state 
Potts model (defined geometrically for Q generic using the loop/cluster expansion) on the 
square lattice. This solution is based on the detailed analysis of the Bethe ansatz equations 
(which involve staggered source terms of the type "real" and "anti-string" ) and on extensive 
numerical diagonalization of transfer matrices. It involves subtle distinctions between the 
loop/cluster version of the model, and the associated RSOS and (twisted) vertex models. The 
essential result is that the twisted vertex model on the transition line has a continuum limit 
described by two bosons, one which is compact and twisted, and the other which is not, with 
a total central charge c = 2 — | , for y/Q = 2 cos j . The non-compact boson contributes a 
continuum component to the spectrum of critical exponents. For Q generic, these properties 
are shared by the Potts model. For Q a Beraha number, i.e., Q — 4 cos 2 ^ with n integer, 
and in particular Q integer, the continuum limit is given by a "truncation" of the two boson 
theory, and coincides essentially with the critical point of parafermions Z n -2- 

Moreover, the vertex model, and, for Q generic, the Potts model, exhibit a first-order 
critical point on the transition line — that is, the antiferromagnetic critical point is not only 
a point where correlations decay algebraically, but is also the locus of level crossings where 
the derivatives of the free energy are discontinuous. In that sense, the thermal exponent of 
the Potts model is generically equal to v = ^. Things are however profoundly different for 
Q a Beraha number. In this case, the antiferromagnetic transition is second order, with the 
thermal exponent determined by the dimension of the ipi parafermion, v — ■^ 2 -. As one 
enters the adjacant "Berker-Kadanoff" phase, the model flows, for t odd, to a minimal model 
of CFT with central charge c = 1— n_i\ t i while for t even it becomes massive. This provides 
a physical realization of a flow conjectured long ago by Fateev and Zamolodchikov in the 
context of Zn integrable perturbations. 

Finally, though the bulk of the paper concentrates on the square-lattice model, we present 
arguments and numerical evidence that the antiferromagnetic transition occurs as well on 
other two-dimensional lattices. 
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1 Introduction 



It is well known how to generalize the definition of the Q-state Potts model to Q real by turning 
to a loop/cluster formulation, and the related six- vertex model (see, e.g., ^IHEl)- On the square 
lattice, the ferromagnetic Potts model thus defined exhibits reasonably well understood critical 
properties. In many respects, it can be considered as the "lattice equivalent" of the Dotsenko- 
Fateev twisted bosonic theory, and thus, has served as a benchmark for our understanding of 
numerous issues in conformal field theory and integrable systems. 

The antiferromagnetic case is much less understood. This is not too surprising, as antiferro- 
magnetic models often exhibit frustration, and thus more complex phase diagrams. While the 
physical interest of this model per se is probably limited to a few special values of Q, it turns out 
that the difficulties encountered in its study are very similar to those encountered in the study 
of models with supergroup symmetries. These models are believed to play a crucial role in the 
description of critical points in non-interacting disordered systems. 

Strong similarities with the s/(2/l) spin chain relevant to the spin quantum Hall effect [I] in 
particular motivated us to tackle the long overdue continuation of a first work on the topic by 
one of us in 1991 2 . In this paper, we are finally able to present a conjecture, supported by 
a long list of arguments and partial calculations, for the conformal field theory describing the 
critical antiferromagnetic Potts model on the square lattice for Q £ [0,4]. In a nutshell, this 
theory involves two bosons: a compact boson not unlike the one appearing in the description 
of the ferromagnetic critical point, plus a non-compact boson. To our knowledge, this is the 
second time only that a non-compact boson is encountered in a statistical mechanics model with 
a finite number of degrees of freedom per site, the first time being the super Goldstone phases of 
jSJ. The appearance of such bosons is generic in the study of conformal field theories based on 
supergroups. 

Another major peculiarity of the antiferromagnetic transition — and of the neighboring "Berker- 
Kadanoff " phase — is the profound role played by the quantum group symmmetry, and the related 
differences between the Potts model and the vertex model for Q generic, as well as the difference 
between the case Q generic and Q a Beraha number [i.e., Q = 4 cos 2 ^ with n integer] for the 
Potts model. 

Finally, the last peculiarity of the antiferromagnetic transition is that it is, for Q generic, a 
first-order critical point, that is a point where correlations are algebraic but at the same time the 
derivatives of the free energy are discontinuous. This feature was initially observed in Ref. [5] for 
the particular case of the Q — > limit. 

Obviously this quick summmary of our conclusions requires much elaboration, and careful 
distinction between incarnations of the model which are usually considered "equivalent". For 
easy reference, we review in section 2 key aspects of the various representations of the Potts 
model. We also comment on the importance of boundary conditions and on the definition of 
the transfer matrices in the various representations. Further background material is given in the 
main text, and details can also be found in [2]. 

The rest of the paper is organized as follows. We shall start in section 3 by the simplest 
problem, which is the properties of the (mostly untwisted) vertex model associated with the 
Potts model on the antiferromagnetic critical line. 

In section 4, we proceed to discuss the properties of the twisted vertex model, as well as those 
of the loop model and thus the Q-state Potts model for Q generic. 

In section 5, we discuss the related incarnation of RSOS models for Q a Beraha number. We 
show how their continuum limit is given by a parafermion theory, and use this in particular 



2 



to complete the identification of the continuum limit of the antiferromagnetic Potts model in 
section 6. 

Another subtlety in this problem is the parity of the number of sites (Potts spins) in the 
transfer matrix. While most of the paper discusses the case N even, interesting things also 
happen at N odd, making the picture richer (the non-compact boson acquires a twist). This is 
discussed in section 7. 

Section 8 goes back to physics, with the discussion of the magnetic exponent and the extremely 
delicate case of the thermal exponent. We discuss there also the issue of the first-order critical 
point, and the flow from parafermionic to minimal theories a la Fateev-Zamolodchikov. 

Section 9 contains conclusions as well as prospects for (much) future work. 

The appendix meanwhile discusses the limit Q — > further, in particular how the R matrix 
for blocks of four vertices in the staggered vertex model case reduces to the osp(2/2) R matrix 
in this limit. 



2 Representations of the Potts model 

The Q-state Potts model is defined initially in terms of integer- valued spins cr, = 1, 2, . . . , Q living 
at the lattice vertices i, and interacting through nearest neighbour coupling, with the partition 
function 

Z = ^[]e M ^. (1) 

Here K is the coupling constant and (ij) denotes the lattice edges. We shall refer to this as the 
Potts spin representation. 

The low-temperature expansion of Eq. Q gives the cluster representation |Hj for which 

Z= Q< E \e K -l)W, (2) 

EQ(ij) 

where c(E) is the number of clusters (connected components) in the edge subset E. The equiv- 
alence between the Potts spin representation and the cluster representation holds true with any 
boundary conditions. Note that now one can decide to give the variable Q any arbitrary real 
value, hence extending the definition of the Potts model to Q non-integer. 

Instead of the clusters, one may count the number of loops 1(E) of the cluster boundaries 
which live on the medial, or surrounding, lattice (here another square lattice oriented diagonally 
with respect to the initial one). Using the Euler relation one obtains the loop representation 

Z = Q S ' 2 £ Q l{E)/2+v{E) [-7^) ( 3 ) 

E c(ij) \ y ) 

where S is the number of lattice vertices. The function r](E) = 0, except on a torus if there is a 
cluster in E that wraps around both periodic directions in the latter case n(E) = 1. With this 
slight topological subtlety, one obtains a loop representation of the Potts model for Q arbitrary, 
totally equivalent to the cluster representation. Formulas (j2J) and (j2J) define the model we want 
to study, and for which we want to determine the phase boundaries, and critical properties. 

A key feature of the Q-state Potts model for Q arbitrary is that it has a non-local definition. 
While this is not a problem for applications (e.g., the study of percolation, trees, forests, etc.), 
it makes the analytical and numerical studies, or the construction of associated conformal field 
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theories, particularly tricky. Let us see this more explicitly. We will be mostly interested in the 
case of the square lattice, which we will often consider as being built up by a transfer matrix acting 
along the "time" (t) direction (of length M lattice spacings), the perpendicular direction being 
referred to as the "space" (x) direction (of width N lattice spacings). We then have S = NM. 

We consider for definiteness the case in which the transfer matrix propagates along 1 the spin 
lattice (i.e., diagonally with respect to the medial lattice), and we let the width of the x-direction 
be N Potts spins 2 . 

The basis states for the transfer matrix in the Potts spin representation, valid for Q integer, 
are of course N variables &i = 1,2, ... ,Q, whence dim T = Q N . In the cluster representation, 
basis states are set partitions of N elements respecting planarity, whence dim T = Cn with 
= n\(n+i)\ being Catalan numbers. For N S> 1, we have CV ~ 4^. Excited sectors are 
constructed by restricting to set partitions with > iVo components, of which exactly Nq are 
marked. By definition, the marked components cannot join among themselves under the action 
byT. 

In the loop representation, basis states are complete pairings of 2N points respecting planarity, 
with again dim T = CV- Excited sectors are constructed by replacing Nq of the pairings with 
exactly 2Nq marked loop segments that cannot close among themselves. 

The crucial point is now that the partition functions (j2J or (jHJ) cannot be written as traces 
of the transfer matrix raised to the power M because of the non-locality. This is of course 
expected, the cluster or loop transfer matrices having different dimension (and thus different sets 
of eigenvalues) than the Potts spin transfer matrix for Q integer. Instead, the partition functions 
are written via more complicated "functionals" of the cluster or loop transfer matrices, alllowing in 
particular for the disappearance of many eigenvalues (corresponding to "non-local" information) 
for Q integer. For instance, with boundary conditions which are free in the x-direction and 
periodic in the t-direction, the loop/cluster partition function is given by Zm = (vf\T M \v\) for 
suitably initial and final vectors (and in fact with a modified transfer matrix defined in terms of 
connectivities among two time-slices) ^0]- The end result is that 

N 

Z M = ^(2iVo + l) g E(4 7Vo) ) M , (4) 

7V =0 k 

where {A^ ^} is the set of eigenvalues of the foregoing transfer matrix, corresponding to the 
sector with Nq marked components defined above. The amplitudes (2 No + l) q are ^-deformed 
numbers and are derived from quantum group considerations [H], while we have set ^fQ = q + q~ l . 

The essential point is that the amplitudes in Q can vanish when Q is integer. In fact, setting 
q = exp(z7r/t), they can vanish whenever t is rational. A remarkable feature of the Potts model 
around the antiferromagnetic transition is that the disappearing eigenvalues can be the leading 
ones, giving rise to strong singularities in the thermodynamic properties. In contrast, in the 
ferromagnetic region, the disappearing eigenvalues correspond only to excited states, and lead 
merely to the disappearance of some non-local observables from the spectrum of excitations. 

To summarize, the lesson so far is that for q a root of unity, not all eigenvalues of the cluster 

1 In some cases it may be advantageous to study instead the case where the transfer matrix T propagates 
diagonally with respect to the spin lattice (i.e., axially with respect to the medial lattice). 

2 In all subsequent figures of numerical data, the finite-size estimates of free energies /, central charges c, and 
physical scaling dimensions x — h + h are labelled by the width N in terms of Potts spins. When the transfer 
direction is diagonal with respect to the spin lattice, we use the label N*, meaning that the projected width is 
\/2N lattice spacings in terms of the spin lattice. 
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or loop transfer matrix do contribute to the Potts partition function. The same would of course 
be true in the doubly free, or doubly periodic case. 

It would now seem crucial to remove the non-local aspect of the weights in the loop/cluster 
representation. To this end, we again consider the case of free boundary conditions in the x- 
direction and periodic in the t-direction. For all loops which are homotopic to a point, assigning 
independently an orientation to each loop and giving a weight q a / 2n to a loop that turns through 
an angle a to the left then reproduces jllj the full loop weight of Q 1//2 upon setting Q 1 / 2 = q+q^ 1 - 
As for the non-contractible loops, they would acquire instead a weight 2 in this formulation, but 
this can be corrected by assigning the weight q (resp. q^ 1 ) to any oriented loop that passes 
through the periodic boundary condition. Summing over the oriented loop splittings of vertices 
which are compatible with given edge orientations finally gives a representation in terms of the 
six- vertex model This model again lives on the medial lattice. The Hamiltonian of the 

corresponding (XXZ) spin chain can be extracted by taking the anisotropic limit, and is useful 
for studying the model with the Bethe Ansatz technique. Note that the vertex model is in general 
staggered, as the weight assigned to a given vertex depends on whether this vertex stands on a 
direct or a dual edge of the original spin lattice. Also, in the case of boundary conditions which 
are free in the x-direction and periodic in the i-direction, the vertex model has special weights for 
the vertices at the surface, corresponding to boundary fields in the XXZ formulation. We shall 
refer to it as a "twisted" vertex model. An untwisted version would be obtained by setting the 
boundary fields to zero. 

Unlike for the cluster or loop transfer matrix, the transfer matrix of the vertex model is a 
local object. 3 Basis states are orientations of 2N arrows, with an equal number of ups and downs, 
whence dim T = {N + l)Cj\r. Excited sectors are constructed by having instead N + k ups and 
N — k downs. Periodic boundary conditions in the x-direction are imposed by giving a special 
weight (twist) to the first arrow in each row. But again, the correspondence with the Potts model 
implies that the object of interest is not the trace of the transfer matrix raised to the power M, 
but some functional of it. In the simplest case of free boundary conditions in the x-direction and 
periodic in the t-direction, adjusting the weight of non-contractible loops leads to the expression 
Zm = tr q s *T M , where S z is the conserved spin (arrow flux) along the transfer direction. The 
commutation of the vertex transfer matrix with the quantum group SU(2) q leads to a more 
compact expression of Zm in terms of g-dimensions, with the same features of disappearing 
eigenvalues when q is a root of unity. So we get the same lesson that for q a root of unity, not 
all eigenvalues of the twisted vertex model transfer matrix do contribute to the Potts partition 
function. The same would of course be true in the doubly free, or doubly periodic case. 

Obviously the phenomenon of disappearing eigenvalues requires great care in handling the 
results from transfer matrix analysis. 

To complete this introduction we must now discuss briefly the doubly periodic case. The 
correspondence between the models is then more complicated. Consider for instance the loop 
representation and its vertex model reformulation. The weig ht Q 1 / 2 for non-contractible loops can 
only be reproduced, after orientation, by giving a special weight that depends on the winding of 
the loop in both periodic directions; this weight, in turn, transforms into a complex combination 
of weighed sectors for the vertex model obtained after summing over loop splittings! On top 
of this, the sector with no non-contractible loops must also be singled out because of the extra 
Q"n( E ) factor in Q ! This is discussed in details in the literature [0] , and will be addressed below. 
For now, we need to consider only what happens in the vertex model when the non-contractible 

3 Note however that the presence of a seam in the twisted version remains a non-local feature. 



5 



loops wind in only one direction. We have already seen what to do when they wind in the t- 
direction: their weight must be adjusted from 2 to Q 1 ^ 2 by introducing a modified trace term. 
Similarly, when they wind in the x-direction, one must introduce a twist term assigning the 
weight q (resp. q^ 1 ) to any oriented loop passing through the x-periodic boundary condition. 
This leads, in the XXZ limit, to a modified coupling between the first and the last SU(2) q spins, 
and defines what we call the twisted vertex model or spin chain transfer matrix. This is mostly 
an intermediate object for us: while we consider the Potts partition function or the untwisted 
vertex model partition function in this paper, we do not enter the discussion of what would be 
the proper definition of the twisted vertex model partition function. Going beyond the simple 
case of loops which are non-contractible in the x-direction is a very technical topic |12| (best 
handled in the framework of Temperley-Lieb and periodic Temperley-Lieb algebras) which we 
will address when necessary in the following. Suffice it to say that for generic q, the eigenvalues 
of the cluster or loop transfer matrix are a subset of the eigenvalues of the vertex model transfer 
matrix for various twist sectors. 



3 Bethe Ansatz analysis of the vertex model transfer matrix 

The antiferromagnetic critical line for the square-lattice Potts model (defined, for Q generic, 
using the loop/cluster formulation) was identified by Baxter |13j as 

(e*i + l) (e K * + l) = 4 - Q (5) 

where K\,Ki are the usual horizontal and vertical couplings. It is as always convenient to 
introduce the variables 

e K - 1 

and set 

sinii 
sm(7 — u) 

where y/Q = 2 cos 7, with 7 G [0, |], and u is a spectral parameter. The antiferromagnetic line 
then corresponds to 

1 cosii smv 7T . , 

— = 7 t = ; n> v = u±- (8) 

X2 cos(7 — u) sin(7 — v) 2 

The Potts model on the line (J5J) is never selfdual. Selfduality would correspond instead to x\ = i, 
that is u = v. 

The Potts model is generally related to a six-vertex model with quantum group symmetry 
obtained from the loop formulation by putting arrows on the medial lattice (see section 2, and 
Ref. [2] for further review), and using the spectral parameters u,v to define the vertex R matrix 
on the two sublattices. In general, this vertex model is not solvable. 

The peculiarity of the antiferromagnetic line is that the vertex model — and hence the Potts 
model, at least in principle — becomes solvable because, when v = it can be considered as a 

particular case of models exhibiting U Z invariance" |14j . In general, such models are obtained by 
associating a spectral parameter hi with every horizontal line (with label i), a spectral parameter 
Vj with every vertical line, such that the R matrix at the vertex ij is R(hi —Vj). On the 
antiferromagnetic line, the Vj are alternatingly 0,^, while the hi are alternatingly u,u + ^. 
Thanks to the special shift §, and the fact that u can be identified with u ± ir, there are indeed 
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only two types of vertices for the resulting vertex model, with interactions given by R(u) and 
R(u + f ) respectively (see the appendix for more details). The isotropic case (K\ = K2, x\ = 22) 
would correspond to n = ^ + |. 

While vertex models with imaginary heterogeneities of the spectral parameter have been 
studied in some detail ^H] , the situation of a real heterogeneity does not seem to have attracted 
much attention. We will now see how the vertex model exhibits truly remarkable properties. 

We consider a transfer matrix propagating in the vertical direction. The basic equations can 
be obtained using standard techniques — they appeared intially in the pioneering paper by Baxter 

m 



and 



3.1 Ground state energy, degeneracies 

We will use modern Quantum Inverse Scattering notations, following [T^]. We introduce the 
kernel 

, , sinh i(a + 217) , , 

0(a) =i In— -21 (9) 

smn 2 (a — ^7) 

together with the source terms 

sinh ^ (a -yy) 

pi(a) = zln , . — — — (10) 
smh^a + i-y) 

, , cosh \ (a — ry) , , 

p- x a =iln f / (11) 
cosn 2 (a + ^7) 

Recall that the parameter 7 is related to Q by y^Q = 2 cos 7, with 7 G [0, §]■ The Bethe equations 
read then symbolically, for a system of 2N vertical lines (this corresponds to N Potts spins), 

f e m(°t)^j N f e iP-l(°i)\ N TJ ^4{ot-af) _ e ii> ^2) 

a' 

Here, we have introduced an adjustable phase ip corresponding to twisted boundary conditions. 

The eigenvalue of the transfer matrix has a complicated expression, but in the limit of N 
large to which we will restrict, it reduces to 

A oc Y[ e* Pl < - Q+2i "^ e iP-i(a+2iu) ^3) 
a 

Rather than the vertex model transfer matrix, we will consider the hamiltonian, obtained in the 
strongly anisotropic limit u — > 0. Its eigenvalues read 

E = cY / -^[p 1 (a)+p. 1 (a)} (14) 

where c is a (positive) constant to be chosen below. 

A remarkable property of these Bethe equations is that they, like the energy, are invariant 
under a global shift of the roots a — > a+iir. This is due to the peculiar form of the heterogeneities; 



4 Note that some misprints have cropped up in this reference: in Eq. (4.40) there should have been a t~ 2 instead 
of t 2 on the right-hand side, as confirmed by the analysis of the analytic Bethe Ansatz, Eq. (5.4) of that same 
paper. 
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in general, the invariance is only under a — > a + 2iir. Equivalently, one can think of this 
as invariance under 7 — * 71 — 7. States for which the root configuration is invariant under 
q — > a + iir will correspond to non-degenerate eigenvalues, while states which are not correspond 
to eigenvalues degenerate twice. 

Writing the contribution of the roots to the energy more explicitly as 

E = cV u 2sln27 (15) 

^ cosh 2a — cos 27 

we can conjecture that the ground state is obtained (at least when the phase ip = 0) by filling 
the lines Im a = \ and Im a = — ^. The source terms are of a new type then, which we will 
denote qi,q-i- The Fourier transform (with conventions f(k) = f(a)e lka da) is 

- - sinti7fc 

«' + ^ = - 2ir imhi^ (16) 

The other kernels to consider are 0n = and 0i,-i(a) = (j){a ± in), with Fourier transforms 

sinh(-7r — 27) k 



sinh7r/c 
sinh 27 A; 



01,-1 = 2tt ' (17) 

sinn ttk 

The minus sign in the source term leads to a ground state where rapidities are a decreasing 
function of the Bethe integers, and we get, in the continuum version 

27r[pi+pi] = + - 011 *Pl - 01,-1 *P-i 

2vr[p_i + p^] = -gi + - 0-i,i *pi - (18) 

Solving these equations gives the densities in the ground state. It is most useful to go one step 
further and write "physical equations" where on the right-hand side appear holes in the ground 
state distributions, that is densities of excitations over the physical vacuum: 

27r(pi+p?) = s(a) + $n* + f£ x 
27r(p_i + p^) = s(a) + $_i,i*^ + $-1-1 *p\ (19) 

The kernels here are defined by their Fourier transforms 

cosh(-7r — 37)^ 



$11 = <3?_i _i = -2ir 



2 sinh jk sinh(7r — 27)^; 



<l $2 cosh 7/0 

2 sinh 7 A; sinh(7r — 27) k 

The function s = — w n r , /n , s(a) = -, — „ t n , — * — . From this the ground state densities first 

cosh(7r— 2~f)k/2 ' \ / (77—27) cosh 7r _ 27 Q 

follow, pf s = p^ = ^rs(a). The ground state energy per vertex is then 

E gs _ _f sinhfrfc) 



N J sinh(7rfc/2)cosh(7r - 27)/c/2 

It is interesting to compare this result with Baxter's expressions for the free energy [T^]. Using 
analyticity techniques, he finds that the free energy per site can be written as / = f{Ki)+f(K2) = 
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f(u) + f(u + vr/2) where the function / (denoted by 4> in his work) is given in Eqs. (30a)-(30d) 
of JH]- To compare with our hamiltonian we need to let u approach from below, and take a 
derivative, so we must consider the expression 



smh 2 ('yk) sinh(2uk) f sinh(7/c) sinh(7r — j)ks'mh.(2uk) 

dnj~ ; — : — ; — ~r — ; — : — ; rr~ ~\~ I dk 



du^^ ) )] ^ u j A; sinh(7rA;) sinh(7r — 27)/;; ' J k sinh(7rA;i) sinh(-7r — 2j)k 

dk sinh( 7 fc) 

sinh(^/2) cosh(^ - 2j)k/2 V ' 

This agrees with our ground state energy, given in (|21|) . and constitutes an interesting check of 
our equations. 

3.2 Excitations 

The "physical" Bethe equations immediately show that there are two branches of massless exci- 
tations (at least) obtained by making holes in the ground state at ±00 on the lines Im a = § and 
Im a = — ^. This is twice the usual result for the XXZ chain, and leads one to expect that the 
central charge will be equal to two, a result we will confirm below. The rapidity of the excitations 
(such that e = ±p tx e 9 ) is 9 = ^jj- a. As for the energy, one has 

N = W + C J S{a) (a) + P - l{a) \ (23) 

so relativistic invariance requires c = — 27 ^ . 

The kernels $ describe a complex scattering theory, which we do not yet entirely understand. 
A crucial fact is that these kernels diverge at small k, a very unusual feature, that is also encoun- 
tered in the studies of models based on super algebras [I]. We will argue below that it implies 
the existence of a continuum component of the spectrum. One can however observe some drastic 
simplifications for "symmmetric excitations", i.e., excitations which are identical for the two lines 
Ima = | and Im a = — f . For these, the "effective kernel" is obtained as $ = <J?n + $1,-1, 

£ sinh(7r - 4'j)k/2 , 

$ = -27r — (24) 

2 sinh 7/c cosh(7r — 2^)k/2 

This is the soliton-soliton kernel for the sine-Gordon model regime with §- = —. This observation 
has a simple origin. Consider again the Bethe equations: it is clear that there is a particular type 
of solutions obtained when the real parts of the roots on Im a = ^ and Im a = — ^ coincide. 
The equations for these real parts are then of the form 

sinh(a — ij + in/2) \ iibTl sinh(a — a' — 2ij) 



while the energy reads 



sinh(a + 27 — iiT/2) ) sinh(a — a' + 2ij) 

E = _^-^y2 2^27 

7r ~ cosh 2a + cos 27 
Comparing with the results for the XXZ chain (with relativistically invariant normalization) 

H = -^jw ? M^ 1 + ~ C ° S 7Vf (2?) 
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we see from standard works on this topic ^H] that this spectrum will lead to twice the central 
charge and critical exponents of the XXZ chain (|27|) with parameter j' = tt — 27. On the other 
hand it is known that the conformal field theory limit of this XXZ chain is a free boson with 
coupling constant g = ^ (in agreement with the value ^- = ^ obtained before for the scattering 
matrix). Notice in particular that the limit Q — > corresponds to 7' = 0, i.e., the XXX chain. 
In other words, there has to be a whole class of excitations in our system for which the gaps 
will be the same as the gaps of this XXZ chain ()27|) . up to a factor two (this should hold with 
or without twist angle ip). We call this 'the XXZ subset' of our spectrum. It is tempting to 
speculate that it corresponds exactly to the configurations of roots which are invariant under the 
extra Z 2 symmetry a — > a + in. Of course, there has to be many more excitations, since doubling 
the spectrum of the XXZ chain creates many gaps in the conformal towers. These excitations 
should appear with an extra degeneracy of two, since they corresponds to configurations of roots 
which are not symmetric. 



3.3 A quick analysis of finite-size effects 

We can make further progress on this question by analyzing in more details the finite-size spectrum 
from the Bethe Ansatz. For this, we go back to the physical Bethe equations involving holes in 
the distributions of strings and antistrings. It is interesting to determine the dressed charges of 
the excitations. Normally, this is an easy task. Write the bare equations as p + p h = p + Kp. It 
follows that p + p h = Zp + (1 — Z)p h , where Z = (1 — K)~ . Standard formulae then exist to 
express the gaps in terms of the matrix Z ^7] : 

h + h = Ara*Z -2 An + d t Z 2 d (28) 
where An is the two-column vector ^ ^ , and Araj the change of solutions of the i th type. 

Similarly, d is the two-column vector ( ^ j , and d{ the number of solutions of the i th type 

backscattered from the left to the right of the Fermi sea. In (|28|) . t denotes transposition. 

Here however we encounter a major snag: the matrix Z has all its elements infinite, as the 
kernels diverge when k — > 0. We have taken the most naive approach possible, and assumed that 
the standard formulas still made sense in this case, after maybe some regularization. It turns out 
that such regularization is easy: just keep k small but non-zero in the matrix elements. One can 
then calculate the inverse of Z, and let k — > in the end. One finds 

*- 1 = »m^ = ?f^ I'l) (29) 



fe^o Z{k) 



and thus 



A + A = -^(Ani + An 2 ) 2 + x (Am - An 2 f + ^-{d x + d 2 f + 00 x (d 1 - d 2 f (30) 

Z7T 87 

We recover here the XXZ subset (see later) for d\ = d 2 . On top of this however, we have a 
continuum of soft modes corresponding to Ani 7^ An 2 . This is similar to observations in ^H] and 
[I] where the soft modes were interpreted as arising from a non-compact boson. We therefore 
believe that the continuum limit of the untwisted vertex model is described by a compact boson 
(called 4>\ in what follows) with g = \ plus a non-compact boson (called 4> 2 ) giving a continuum 
of soft modes on top of the gaps for the first boson. 
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Figure 1: A block of four vertices can be considered as a "supervertex" for a model with 2x2 = 4 
degrees of freedom per edge. 

While the action of the twist angle tp on <pi (that is, the XXZ subset of the spectrum) is well 
understood, the action on fa seems to be more subtle, and cannot be captured by our simple 
analysis. We shall see that even the parity of N affects 4>2 hi a non-trivial way. 

3.4 Some remarks on extra symmetries 

The Bethe equations expressed in terms of roots at Im a = ±^ look like equations for a rank two 
algebra, and suggest the existence of two £7(1) conserved quantities, that is, an extra quantity on 
top of the spin S z . One can get an idea of what this extra quantity might be by considering the 
limit Q — > 0. We observed in JH] that in this limit, the Bethe equations and the energy of the 
model coincide with those of the integrable osp(2/2) spin chain in the fundamental representation. 
One can in fact go further (see the appendix), and prove that the osp(2/2) integrable R matrix 
can be obtained by considering a "block" made up of four elementary vertices of the 6 vertex 
model. In this case, the extra U(l) charge is the B quantum number of the osp(2/2) algebra; 

this is defined initially for a pair of neighbouring spins \S2i-iS2i) by B = for | + +) or | ), 

and B = ±1/2 for | H — ) ± i\ — h) (see Fig. and then extended to the whole system by 
summing over i. It is very interesting to study what becomes of this extra symmetry on the 
whole antiferromagnetic line, and the relation with the osp q {2/2) integrable R matrix — we will 
report about this elsewhere. In any case, the quantity Ani — An2 plays the role of an extra U(l) 
charge along the whole antiferromagnetic line. 

On top of this charge, the antiferromagnetic line exhibits an extra Z2 symmetry, as already 
mentioned in 2 . This can be seen most clearly by considering again the block of spins in Fig. ^ 
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and writing its R matrix algebraically as (for notations see [2]) 

R = X 2 (u)X 1 (v)X 3 (v)X 2 (u), « = u±| (31) 
Using the Yang-Baxter equation repeatedly, one can show that the operator 

commutes with R, where the fact that X (^) = X (— ^) is crucial. On the other hand, one has 

X (±0 = P x - P (33) 

where Pi,Po are the projectors on (quantum) spin one and zero representations in the product 
of two spin 1/2 representations, so C 2 = 1. 

We can now discuss briefly the phase diagram of the Potts model in terms of the symmetries 
of the associated vertex model. The selfdual case x\ = gives rise to two kinds of physical 
behaviors, which correspond in the isotropic case to two selfdual lines, x% = x% = 1 and x\ = 
x 2 = — 1. In the selfdual case, the symmetry of the associated vertex model is £7(1) X Z 2 . The 
£7(1) is just the conservation of the spin S z , while the Z 2 is the symmetry of translation by 
one lattice site. The latter symmetry, in the continuum limit, becomes another £7(1) JH] so the 
symmetry of the continuum limit on the selfdual lines is £7(1) x £7(1): this corresponds to chiral 
symmetries for left and right fermions, or symmmetry (pi — > <pL + a, <j}R — * + (3, of the free 
bosonic continuum limits. Going away from the selfdual lines breaks the Z%. If one starts from 
the physical critical line (x\ = x 2 = 1), the model becomes massive. If one starts from the 
non-physical critical line {x\ = x 2 = — 1), one enters the Berker-Kadanoff phase [2] where the Z 2 
symmetry breaking operator is irrelevant, and the full £7(1) is actually restored in the continuum 
limit. 

The antiferromagnetic critical line marks the boundary of the Berker-Kadanoff phase. As 
just discussed, on this line an extra £7(1) x Z 2 symmetry appears, and it is reasonable to expect 
that in the continuum limit, this also gets enhanced into the full £7(1) x £7(1), corresponding 
to another free boson. This is in agreement with the Bethe Ansatz calculation of the central 
charge. However, when Q = it has been argued in ^Hj that this boson is non-compact. This 
feature must be true along the full antiferromagnetic line in fact, as follows from the finite-size 
corrections to the gaps. 



3.5 The spectrum of the untwisted vertex model for /V even 

One subtlety is that results for N odd and N even turn out in fact to be profoundly different. We 
do not understand this fully, but notice this can be expected, as the foregoing analysis applied 
mostly to N even, which allowed one to treat the two lines Ima = The case of N odd will 
be considered later. In what follows, we parametrize 7 = j, with t E [2, oo]. 
To start, we consider the XXZ subset some more. Magnetic gaps read 

A = 2xS Jtl = (i!l (34) 

where S z is the spin of the "equivalent XXZ chain" — that is the chain for which the Bethe 
equations are (|25|). Since the antiferromagnetic system has 2N spins, whereas the XXZ chain 
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described by (f25|) is for N spins only, we must rescale the spin according to S z = where S z 
is the spin in our initial system, i.e., S z = Arti + An2 in 1|30|) . Therefore 

A = ^2 (35) 

where S z (the spin in our initial system) is an integer since we have a even number of spins. 

One can see here again why parity effects are likely to occur. If N is even, since we must 
also have an even number of upturned spins to make up pairs of Bethe roots in the "folding" of 
the equations, the result strictly speaking will only apply to S z even, including the ground state. 
Numerical study shows that in fact it holds for S z odd as well. If N is odd on the other hand, 
since we must also have an even number of upturned spins to make up pairs of Bethe roots in the 
"folding" of the equations, the result strictly speaking will only apply to S z odd. Numerically, 
we have observed indeed that for S z even, including the ground state S z = 0, results for N odd 
obey a different pattern. See section 7 below. 

The electric charges in the S z = sector of the equivalent XXZ chain are integer (e = d\ = c^), 
and give rise to the following gaps A = 2 x |^ = Combining all the excitations, we thus get 

exponents from a single free boson, A em = | ± e\/t\ , with e and m integers. 5 

Meanwhile, the non-compact boson adds one unit to the central charge, and 'decorates' these 
gaps by soft modes, which would manifest theshelves in a numerical study as very large de- 
generacies One can write this in a compact form by considering the generating function of 
levels 

Zcon] = £ g A-c/24 -A-c/24 = f fo^-l/M x ]_ £ q A em -A em ^ 
A, A ^ e > m 

with A em (A em ) = | (jj^ ± e*Jt\ . It would coincide with the partition function of the doubly 

periodic vertex model with q = e - 2nM / N the modular parameter for a system of size N x M and 
rj = gV24 n£° =1 (l - q n ) Dedekind's eta function. 

This means that in the sector of vanishing magnetization, one should observe a continuum 
spetrum above the ground state. In the sector with magnetization unity, one should observe a 
continuous spectrum over the basic magnetic gap, etc. Numerical study by direct diagonalization 
of the transfer matrices (possible for N < 12) is compatible with this result, though the sizes 
accessible are too small to provide a complete confirmation. In [3] where the Bethe equations 
themshelves were studied for a simimlar problem, it was found that sizes as large as N ~ 10 4 
were necessary. 



4 The spectrum of the twisted vertex model and the Potts model 
transfer matrices on the antiferromagnetic line 

The usual strategy is, once the continuum limit of the untwisted vertex model is identified, to 
evaluate the effect of the twisting and other ingredients involved in the correspondence with the 
loop/cluster formulation, and gain in this way knowledge of the conformal properties of the Potts 
model itself. The process in the antiferromagnetic region is however more involved than in the 
ferromagnetic region. To start, we discuss the Berker-Kadanoff phase in more details. 

5 In general, we will denote by A the (rescaled) gaps with respect to the central charge c = 2, and reserve the 
notation h for the (rescaled) gaps with respect to the true ground state of the antiferromagnetic Potts model. 
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4.1 The Berker-Kadanoff phase 



Recall first that the untwisted vertex model in the whole Berker-Kadanoff (BK) phase is described 
in the continuum limit, like on the non-physical selfdual critical line, by a single free boson with 
g = j. It is important to notice that this free boson is in fact exactly the same as the one 
describing the XXZ subset on the antiferromagnetic critical line. 

The first step in going from the vertex model to the loop/cluster representation of the Potts 
model is to introduce a twist in the vertex model. We have explained that this twist has a con- 
venient geometrical interpretation if one decomposes the vertex configurations in loops: without 
it, loops which are non-contractible around the cylinder get a weight 2, instead of the required 
y/Q for the loop representation of the Potts model. Suppose we want to give to non-contractible 
loops a more general weight, w = 2cos7T(X We introduce a seam or twist in the vertex model 
and the associated conformal weight should give the effective central charge 



Since non-contractible loops come in pairs, we have a defined modulo an integer: these are 
the allowed values of the electric charge in the Coulomb gas. The Q-state Potts model itself 
corresponds to a = eo = 7, leading, in principle, to the central charge c = 1 — |. 

This latter value is definitely the central charge of the twisted vertex model. But it is not the 
central charge of the loop/cluster Potts model! 

The point is that the eigenvalues of the loop or cluster model are in general (as discussed 
in section 2) a subset of the eigenvalues of the vertex model, essentially those associated with 
"highest weights" only. The determination of this subset is based on algebraic considerations, 
and follows some simple rules (Hj. For instance, eigenvalues of the vertex model transfer matrix 
in the sector of spin S z = 1 and with no twist are also exactly observed in the sector S z = and 
with a twist a = \ (i.e., giving to non-contractible loops the weight \/Q, that is a = eo). Such a 
"descendent" eigenvalue is thrown away in the loop or cluster formulation. 

In the BK phase, the dimension associated with electric charge eo = 7, since it is the same as 
the dimension coming from the sector with no electric charge (twist) and S z = 1, and equal to 
A eQj o = z?> is thrown away in the truncation from vertex to loop (or cluster). This means that 
the central charge of the twisted vertex model and the loop/cluster model are different in this 
phase. For the vertex model one has c = 1 — f while the central charge of the loop or cluster 

model follows from the next choice, eo = 1 — j, leading instead to central charge c = 1 — 
all this for t generic: this is illustrated in Fig. |2j 

We have performed careful numerical checks of this phenomenon, with some results presented 
in Figs. 13121 A technical note here: one can essentially write two transfer matrices, one propa- 
gating along the axial direction of the original Potts lattice (and thus diagonally for the medial 
lattice, which is the lattice where the loops or the vertex spins live), and one propagating diag- 
onally for the Potts lattice (and thus axially for the loops and vertex spins). Results will not 
(except for parity effects) depend on the geometry in the thermodynamic limit, but might con- 
verge better in finite size with one choice or the other. Note that the first formulation is better 
suited for using the Temperley-Lieb algebra and quantum group symmetries. It is the one we use 
most. In all figure captions we will use the name diagonal and axial in reference to the medial 
lattice (see also footnotes ^£3) • 

Notice that for t an integer (and only then), the eigenvalue associated with eo = 1 — \ is also 
thrown away because it is a descendent from the sector with no twist and spin S z = t — 1, as 




(37) 
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Figure 2: A schematic representation of the behavior of some levels (in finite size) of the vertex 
model in the BK phase, for t generic. Levels which are at the head of an arrow are thrown away 
in going from the vertex to the loop or cluster formulation. The level A = i is thus always 
discarded, leaving generically the next level A = (t — l) 2 /4i to determine the central charge 
c = 1 — 6(i — l) 2 jt for the loop or cluster model for t generic. 




Figure 3: Numerical estimates of the central charge for the twisted vertex model on the non- 
physical selfdual line. The continuous curve is the expected value, c = 1 — 6/t. Diagonal geometry. 
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Figure 4: Numerical estimates of the central charge for the loop model on the non-physical 
selfdual line. The continuous curve is the expected value, c = 1 — 6(i — l) 2 /t. Axial geometry. 

illlustrated in Fig.|SJ This is the mechanism that, in the periodic case, produces the singularities 
of the BK phase at the Beraha numbers Q = 4 cos 2 j, with t integer. It is the equivalent of the 
slightly simpler mechanism studied in [2] in the case of open boundary conditions. In fact, many 
more levels are thrown away, so many that even the free energy per unit site of the (twisted) 
vertex and loop or cluster models are different at the Beraha values in the BK phase. The 
difference of these free energies vanishes as the antiferromagnetic line is approached. Right on 
the antiferromagnetic line, the free energy does not show any singularity while one crosses a 
Beraha number (see section 8, in particular Figs. I2T1422|) , 

A note of caution on the case t integer is in order. Of course, the spectrum of the loop/cluster 
transfer matrix does not show any particular behavior when one crosses these values. What 
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Figure 5: A schematic representation of the behavior of some levels (in finite size) of the vertex 
model in the BK phase, for t an integer. Now the two levels indicated on the right are both 
discarded. Many more levels are in fact discarded, leading to the singularities at the Beraha 
numbers within the BK phase. 
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Figure 6: Numerical estimates of the central charge for the twisted vertex model on the anti- 
ferromagnetic critical line. The continuous curve is the expected value, c = 2 — 6/t. Diagonal 
geometry. 

happens however is that the partition function of the model on a torus, say, would be expressed 
by weighing these eigenvalues with complicated non-integer 'degeneracy factors', and some of 
these factors can vanish when t is integer, a feature we refer to as a level being truncated, or 
discarded. Hence what one calls an exponent of the loop/vertex model has to be defined with 
great care! 

4.2 Central charge on the antiferromagnetic critical line 

We can now carry out a similar analysis for the model on the antiferromagnetic critical line. To 
give non-contractible loops a weight w = 2 cos tree, we introduce a seam, and the effective central 
charge follows immediately as 



leading to a central charge c = 2 — ~. 

Remarkably, this is not only the central charge of the twisted vertex model, but aiso the 
central charge of the loop/cluster Potts model on this line, as was established earlier and can 
be checked numerically to great accuracy. Some results are presented in Figs. EHZI For Fig. El the 
transfer matrix propagates along the diagonal of the medial lattice, for Fig. [7| along the medial 
lattice itself. Similar results are obtained for other geometries. 

This result is remarkable, because we have seen that the central charge of the twisted vertex 
model comes entirely from the XXZ sector, and that this sector is described by the same free 
boson as in the Berker-Kadanoff phase. Hence one should have expected the level corresponding 




(38) 



since now g = — = | as in section 3.2. For the Q-state Potts model itself we have a = bq = j 



t ' 
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Figure 7: Numerical estimates of the central charge for the loop model on the antiferromagnetic 
critical line. The continuous curve is the expected value, c = 2 — 6/t. Axial geometry. 

to (|38|) to be truncated when going from the vertex to the loop/cluster formulation, and the 
central charge of the Potts model to be given by another value. What happens is the following. 

Observe that the XXZ subset of the spectrum (in the scaling limit) on the antiferromagnetic 
critical line coincides with the spectrum (in the scaling limit) in the BK phase, with or without 
twist angle (and up to an additional contribution of one to the central charge). Since, because of 
the doubling effect discussed previously, the vertex model on the antiferromagnetic critical line 
has more scaling gaps than in the BK phase, it means that some of the levels which were at 
finite distance from the ground state in the BK phase have to cross and become scaling levels 
on the antiferromagnetic critical line. (A similar crossing phenomenon is clearly visible on the 
finite-size levels of the loop/cluster model, see Fig. |HJ). One of these levels happens to merge 
(asymptotically) with the descendent level, so after truncation one of them remains, to give the 
central charge c = 2 — 6/t still. This is illustrated schematically in Fig. |5J 

This is confirmed numerically: we find that the watermelon operator h,2 = 0, that is the 
largest eigenvalue of the loop transfer matrix in the sector with zero legs, and weight y/Q for 
non-contractible loops is asymptotically degenerate 6 with the largest eigenvalue in the sector with 
two legs, and weight 2 for non-contractible loops. Since the latter eigenvalue is also observed (as 
a "descendent eigenvalue" ) in the vertex model in the sector with S z = and twist angle a = 1/t 
(giving to non-contractible loops the weight \fQ), it follows that the latter sector has a ground 
state degenerate at least twice. This is the hall mark of a first-order phase transition, and we 
conclude that the entire antiferromagnetic critical line is in fact a line of first-order critical points 
for the twisted vertex model. This is also true of the loop model, even though it has fewer 
eigenvalues, because of the result that /12 = 0. 7 Finally, it is also true for the untwisted vertex 

6 More precisely, the logarithm of their ratio vanishes faster than 1/N. 

7 The existence of a first-order critical point in the loop/cluster formulation was made initially for the Q — > 
limit in Ref. 0- 
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Figure 8: First scaling levels, normalized as free energy densities, of the cluster model with 
t = 5 and size N = 6 (diagonal geometry). The abscissa is the temperature parameter x = 
Q^ 1 ^ 2 (e K — 1). The vertical line shows the bulk antiferromagnetic transition temperature. The 
ground state on one side of the transition becomes a high excitation on the other side, and vice 
versa. 



S z = n 



A=l?L 
At 



S z = 



2 

A=?f- 
4t 



a = 



a = 



Figure 9: The levels corresponding to Fig.|21 but on the antiferromagnetic line. The ground state 
in the sector S z = 0, a = ? has now double degeneracy, so that one of the two levels remains 
after the descendent levels are discarded, leading to the central charge c = 2 — |, for any t. The 
additional level comes down from the very excited part of the spectrum as one moves within the 
BK phase towards the antiferromagnetic critical line. 
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Figure 10: Spectrum of the loop/cluster Potts model transfer matrix in the sector with no 
marked loops/clusters. Two exponents are clearly identified (solid lines) as well as descendent 
levels (dotted lines). Axial geometry. 



model (since its ground state is in fact infinitely degenerate due to the continuous component). 

4.3 Spectrum of the loop/cluster model transfer matrix with no marked loops 

More generally, the weight w of non-contractible loops is compatible with charges e = ± 7 ± n, 
with n an integer. The choice e = eo — 1 = j — 1 leads to a conformal weight (defined over the 
ground state with c = 2 — ~) 

(i-l) 2 -l t-2 

h = ~ 7 = —7- 39 

At 4 y ' 

This is well observed numerically as seen in Fig. I1UI We believe that there is a continuum 
spectrum on top of this gap, all of which is truncated when one goes from the Potts transfer 
matrix to the Potts partition function when t is integer. 

One can also identify numerically a level very well matched by the formula h = |5§ • This has 
no natural explanation within the XXZ subset, and will be discussed in details below. 

We will argue later that the spectrum is in fact continuous above h = This is all for t 
generic of course. For Beraha numbers, the continuous component disappears entirely, as well as 
other exponents. 



4.4 Spectrum of the twisted vertex model transfer matrix in the sector S z = 

Exponents are shown in Fig. 1111 Observe how one sees considerably more levels than in the 
loop/cluster model case (see Fig. I10|) . We will argue later that here, the spectrum is in fact 
continuous above the ground state h = 0. 
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Figure 11: Spectrum of twisted vertex model transfer matrix in the sector S z = 0. Axial geometry. 
Note the considerable set of additional levels when compared to Fig. I1UI The particular levels 
identified with solid lines in Fig. ^] are also present here. 



4.5 Watermelon operators 

These operators describe the properties of L marked lines in the loop model (or, for even L > 2, 
equivalently by L/2 marked clusters in the cluster model). As usual, their conformal weight is 
obtained by considering the ground state of the sector with S z = ^ (L even), so, with respect to 
the ground state of the theory with central charge c = 2 — | , we find 

u, = ^=1 m 

This is also well checked numerically, as can be seen in Fig. 1121 
Notice in particular that the case L = 2 corresponds to /12 = 0. 



5 RSOS versions at the Beraha numbers 

As commented earlier, the Potts model does present singular behaviour at the Beraha numbers 
in the Berker-Kadanoff phase because the largest eigenvalues of the loop transfer matrix do 
not contribute to the partition function. Eigenvalues which actually do contribute are related, 
though not identical, to those of another version of the model, of the restricted solid-on-solid 
(RSOS) type. These RSOS models are particularly interesting to study as they are completely 8 
local models of statistical mechanics, and thus should correspond to well defined conformal field 
theories, presumably "minimal" with respect to some chiral algebra, and hence easier to identify. 

RSOS models can be defined whenever t (in q = e J7f//t ) is rational, t = r/s. Their variables 
are integer heights hi = 1, 2, . . . , r — 1 living on the union of direct and dual lattice sites (i.e., 

8 I.e., they do not even possess a twist. 
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Figure 12: Watermelon exponents on the antiferromagnetic critical line, in the loop model with 
L = 2,4, 6, 8 marked lines. Axial geometry. 

again a diagonally oriented square lattice, but shifted with respect to the medial lattice), subject 
to the RSOS constraint \hi — hj\ = 1 for neighboring i and j. The Boltzmann weights are most 
easily defined by writing the row-to-row transfer matrix |2L)l [2] 

N ( K _ i \ N / K i \ 

* = Q m n + e*-x) n (l + -qIJ^ J (41) 

in terms of Temperley-Lieb generators e,-. The latter are represented as 

(S S i) l l 2 

{ej) hjh , = SQij-u h j+1 ) J] S(fH, Hi K> , (42) 

where S*^ = sin(7r hs/r). 

In the transfer matrix study of the RSOS models, the basis states are collections of 2iV heights 
h±, h±, ... , h2N belonging alternately to the direct and dual lattice, and subject to the constraint 
\hi — hi + i\ = 1. To have periodic boundary conditions in the x-direction we impose h2N+i = hi. 
One can show that dim T ~ [2cos(7r/r)] 2Ar = Q N when N S> 1. 

Of course, the correspondence between the partition function of the RSOS model, the partition 
function of the Potts model, and the spectrum of the loop/cluster transfer matrix or twisted 
vertex model transfer matrix is rich and intricate. For instance, the partition function of the 
RSOS model with doubly periodic boundary conditions on a torus would read f2jQ] 




Z = Q s / 2 £ Q V(E)-l (m/2 ( ) E ( 2 cos ^ ) (43) 
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where the number of non-contractible loops Iq(E) has been singled out. This formula indicates 
that, in the vertex model, r — 1 sectors with different values of the twist are now necessary. 

In this paper we shall restrict the study to t integer (i.e., r = t and s = 1), and use the 
identification of the associated conformal field theory to strenghten our general results (on the 
critical ferromagnetic line for instance, the CFT for the RSOS models would be the minimal 
theory with central charge c = 1 — Q/t(t — 1)). 

The critical antiferromagnetic variety does not intersect the selfdual manifold, and thus cor- 
responds to a staggered RSOS model. We are not aware of any previous study of such model. 

To proceed, let us define the four regimes 
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The Potts model on the antiferromagnetic critical line involves two kinds of vertices, with spectral 
parameters differing by ^. Hence, either they lie in regimes (i), (Hi) or in regimes (ii), (iv). 

We can go from the Potts model to the RSOS model algebraically by using a different rep- 
resentation of the Temperley-Lieb algebra (that is, 'quantum group restricting' the vertex model 
representation). A homogeneous spectral parameter u (such as would be obtained on the self- 
dual lines) would lead to a homogeneous RSOS model of the type studied by Andrews, Baxter 
and Forrester T° dispel confusion, we will reserve the name ABF to homogeneous models, 
and denote the restricted version of the vertex model on the antiferromagnetic line (as well as 
elsewhere in the phase diagram) as the RSOS model (in general, staggered). Where, in the ABF 
phase diagram, these models stand is an interesting question. For u in regimes (i) or (Hi), the 
ABF model is at the transition point between its regime III and regime TV (capital letters refer 
to the regimes originally defined in |U), while for u in regimes (ii), (iv), the ABF model stands 
at the transition point between its regime / and regime II. Based on this observation, we see 
that we ought to expect that the antiferromagnetic Potts model, for a given value of Q, can be in 
two different universality classes, depending on the anisotropy. The isotropic case is what we are 
interested in here: K\ = K2 corresponds to the SDP (self-dual Potts) models in regimes (ii), (iv) 
and thus the ABF models at the transition between regime / and TT. In fact, only the regime 
(ii) corresponds to a physical ABF model, and thus was studied in |21j. 9 

Now remarkably, the central charge we have found for the isotropic Potts model on the 
antiferromagnetic critical line coincides with the one of the associated ABF model at the transition 
between regime I and regime II for t an integer, t = k + 2. It is therefore tempting to conjecture 
that the RSOS version of the antiferromagnetic critical line for k = t — 2 an integer is in the 
universality class of the ABF models at the transition between regime I and regime II, and thus, 
following 121 [ \22\ is a theory of parafermions. Let us discuss further the evidence for this. 

An immediate objection to this claim might be that the ABF models are not at a first-order 
critical point. Indeed, they are not, and the point is that when one goes from the twisted vertex 

9 Note that we also predict from this discussion the existence of more critical regimes for the selfdual Potts 
model itself. Since only the isotropic case was usually considered, attention focussed on the cases x\ = x% = 1 
and xi — xi = — 1, corresponding respectively to u = ^ and u = ^ + f , points which lie in regimes (i) and (iii) 
respectively. So the excursions of the selfdual manifold away from the isotropic points, in regimes (ii) and (iv) 
respectively, remains to be studied. We suspect that the corresponding universality classes will be closely related 
to the ones of the critical antiferromagnetic (hence non-selfdual) Potts model. 
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Figure 13: Same as Fig. |S]but for the RSOS model. Remark the absence of levels responsible for 
the BK phase. Some levels are identical with those of the cluster model, including the ground 
state on the high-rc side of the transition. This ground state does not participate any level 
crossings. 



model to the RSOS models, the eigenvalues from the untwisted S z = 1 sector are discarded 
(L = 2 watermelon), as well as the descendent eigenvalues in the S z = sector with twist a = 4. 
Therefore, the second h = h,2 = exponent disappears, and there is no reason to expect a 
first-order phase transition any longer. This is illustrated in Fig. 1131 

Note that in the ABF language, the extra Z<i symmetry which appears on the antiferromag- 
netic critical line corresponds to the chiral symmetry of the models discussed in |22| . that 
exchanges clockwise and counterclockwise. 

The ABF models at the transition between regime I and regime II have parafermionic expo- 
nents, of the general form [23] 



1(1 + 2) m 2 

h = -j- — ^--rr, l = 0,l,...,k; -l<m<l; I - m = mod 2 



(45) 



and 



(k-l)(k-l + 2) (k-mf 



4(fc + 2) 



4A: 



I = 0, 1, . . . , k; l<m< 2k-l-2; l-m = mod 2 (46) 



All these exponents (and only these) are indeed observed numerically in the RSOS model on the 
antiferromagnetic critical line: see Fig. 1141 

On the other hand, by the quantum group construction [Hj, the eigenvalues of the RSOS 
models are obtained by combining sectors of the vertex model with vanishing spin and twists 
a = jj~^2) ( recau ) t = k + 2), I integer, < I < k. The associated dimension in the XXZ part 
of the spectrum is obtained by setting e = jjQ, and coincides with the one from the sector with 
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Figure 14: Exponents for the RSOS model on the antiferromagnetic critical line, and comparison 
with the spectrum of the model. Solid lines show the exponents (|45|) . while dashed lines 
represent descendents. Diagonal geometry. 

vanishing twist and S z = I + 1; as such, this eigenvalue disappears from the loop model in the 
BK case. On the antiferromagnetic line however, this value is still observed: like in the case 
I = 0, there is double degeneracy in the vertex model, so the loop model still has this exponent. 
Observe now that this exponent (with respect to c = 2 — |) agrees with formula (j45j) for the same 
I and m = 0. 

For a given I, the parafermionic tower gives wei ghts (gSJ). We recognize there the value coming 

2 

from the XXZ subset, from which the quantity ^ has been subtracted. In particular, in the 
sector with twist a = the lowest lying excitation should not be given by m = (unless 
I = 0) but by the dimension of the order parameter in the Zk theory, which corresponds to I = m 
and reads 

»• = (47) 
2k(k + 2) y ' 

Replacing by the value of a gives the result 

((fc + 2)a-l)(fc + l-(fc + 2)a) 

ha ~ 2HkT^) (48) 

The simplest way to put these results together is that the spectrum should contain, for any value 
of a, the XXZ value (with respect to the ground state with c = 2 — |) h = . But on 

top of this, it should also contain other eigenvalues which have double degeneracy, (|48|) being the 
lowest one for a = tt^. Only for / = 0, i.e., a = ^2 wm the two coincide. This pattern should 
presumably extend to t generic. 

Numerical study for the leading exponent in the sector of twist a as a function of a is 
illustrated in Fig. El One sees that for \a\ < h, the XXZ formula holds, while beyond this, it 
has to be replaced by the formula (|48|) with t = k + 2, and t now arbitrary, as expected. 



25 




Figure 15: Effective central charge for t = 6 in the sector with S z = and twist a. For a < j, 
the XXZ gap determines the effective central charge c e g = c — 24/i = 2 — Qta 2 (blue line) , while 
beyond this value, the effective central charge is determined by a different gap, whose origin will 
be related to a non-compact boson later, c e fj = c — 24h a = — 1 + -^{2a — l) 2 (cyan line). Axial 
geometry. 
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Exponents with —l<m<l should also be present, but we have had difficulties identifying 
them numerically in general. An exponent which we have identified however corresponds to (|46[) 
with I = and m = 2: it is the exponent of the first (most relevant) parafermionic chiral current, 
with 

h = ±=1 (49) 

While we have argued that the spectrum has a continuous component for the vertex model 
and the loop/cluster model, there is of course no such component in models, and therefore all 
the corresponding exponents must disappear in the RSOS partition function. This is in particular 
the case for the continuous branch in the sector I = 0, which starts at h = = -A-k. 



6 A free field representation for the antiferromagnetic Potts 
model 



To put all these elements together in the construction of a fully consistent conformal field theory 
for the cluster /loop Potts model does not seem straightforward, and is plagued with the usual 
difficulties, in particular logarithmic features — such a fully consistent construction is not available 
even in the much simpler case of the ferromagnetic critical line. Like in the latter case, some 
progress can be made however by using a free field representation. 

This free field representation is inspired by two features. First, we have deduced the presence 
of a compact boson and a non-compact boson from the Bethe ansatz analysis. Second, the value 
of the central charge at the Beraha numbers t = k + 2 coincides (as was observed in [2]) with 
the central charge for models, and for the latter, a free field representation based on a pair of 
bosons, one compact and one non-compact, is well known [23j (see also |25( [26 ) . 

We thus start, following by introducing a pair of bosonic fields 4>\ and 4>2, with propagators 

(Mz)Mw)) = -2\n(z-w) 
(Mz)Mw)) = -2ln(z-w) 

(50) 

and a stress tensor 

n2 , .• .,2. 



(51) 



With a charge at infinity a® = for the first boson cf>i, the central charge is c 

There are various ways to introduce screening operators in this theory. The first choice, which 



2-1- 



was used in 21 uses both bosons 



J 



2, and leads to the three currents 
l = dfo exp [2ia (pi] 



together with 



J± = exp 



--VUi ± -Vt 

2 2 



202 



(52) 



(53) 



These screening operators, having conformal dimension one, commute with the Virasoro algebra. 
It turns out that they also commute with the "parafermionic" currents 




i + id<p2 exp 



i - iO(p2 exp 



1 



1 



(54) 
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The vertex operators 

Vim = exp 



2 Y / t" 1 ' 2^2" z (55) 
then play a special role: the action of J± on V[ m is only well defined for I integer and I ± m even. 
Vi m is annihilated by the corresponding charges Q± iff — Z < m < I. Similarly, if we consider 
the action of powers of Q±, it is well defined only in the case Q 1 ^ 1 acting on V/ m , unless t is 
rational — we will suppose it is not the case here. Then Q[ +1 V tm = 0. Finally, one can get outside 
the range — / < m < I by acting with the parafermionic fields, which are also annihilated by Q± 
and Q\. 

The case I = leaves fields obtained by successive action of the parafermionic fields on the 
identity, the lowest of which has dimension h = |5§ and is twice degenerate. 

We call the fields obtained in this way the first Coulomb gas contribution. In this sector, (f>2 
presumably does not contribute a continuous part, as the corresponding charges are constrained 
by the screening operators. 

Note that we can also consider twisting differently the boson (f>i with do = |^=. Introducing 
the usual screening charges for a one boson theory, a + a__ = — 1, a + +a_ = 2cto, the corresponding 
screening currents read 

J± = exp [ia±^i] (56) 



The usual Felder construction [27] gives the allowed vertex operators 

$ rs = exp 



r . 1 — s 
— QL+ + i a_ 



(57) 



The central charge and the dimensions are not those expected; c = 2 — 6 ^ p (the additional 

contribution of one coming from the boson 02 ) and h = ^ r ~ s ^ 4 ~^~ 1 ^ . But this theory is non- 
unitary, and admits in particular the negative-dimension operator with r = 0, s = 1, leading to 
an effective central charge c = 2 — I which coincides with the central charge of the loop/cluster 
Potts model. The exponents with respect to this effective central charge are 

h = (tr " g)2 - 1 (58) 

and in particular, the choice r = s = 1 gives the gap h = ^-f L - The -L-leg operators correspond 
meanwhile to /io,L/2 {L even in the Potts model). 

In this sector, cp2 is not constrained, and presumably contributes a continuous spectrum, 
starting therefore right at the effective central charge. 

The "full Coulomb gas" is expected to describe the twisted vertex model: in particular, the 
degeneracy two of the effective ground state with c = 2 — | is reproduced in our picture. Some 
of the levels of the twisted vertex model were represented earlier. 

Going to the loop/cluster Potts model involves discarding some levels, such as the r = 0, 
s = 1 field above. The continuum then should only start at the true central charge, i.e., the gap 
h = with respect to c = 2 — |. 

When t is rational, going to the RSOS model involves discarding many more levels, and 
restricting these Coulomb gases to minimal conformal field theories. In the special case t = k+2 an 
integer, the corresponding minimal model 10 for the second Coulomb gas is M/% + 2,i x [non-compact 
boson <j)2\ and is actually empty. This leaves only the first Coulomb gas, whose restriction gives 

10 We denote by M PtP i the minimal model of CFT with central charge c = 1 — £i£^£J_. 
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Figure 16: The central charge of the untwisted vertex model on the antiferromagnetic critical 
line for N odd takes the value c e g = \ in contrast with the value c e g = 2 for N even. Diagonal 
geometry. 

the Zk parafermionic theory. This agrees in particular with our prediction that the continuous 
component should disappear at t = k + 2 an integer. The case t rational meanwhile should retain 
some of the continuous component. 

7 The case N odd 

Numerical study shows that the ground state in the sector with N odd does not scale like the 
ground state for N even. 11 Maybe the best way to understand this is to start by discussing the 
untwisted vertex model, and the associated watermelon operators. We find that the ground state 
of the untwisted vertex model exhibits an effective central charge 

c cS = ~ (59) 

independently ot the value of t, as illustrated in Fig. (|TH|) . This is in contrast with the value of 
c = 2 we found earlier for the case iV even. 

The simplest interpretation of this result is that the odd parity acts as a disorder line, corre- 
sponding to a disorder operator with dimension A = j^. This in turn can be interpreted within 
the Coulomb gas by assuming that one of the two bosons sees antiperiodic boundary conditions 
in the space direction. The determination of the watermelon exponents (see below) gives the 
values 

(L/2) 2 -l 1 , . 

h L = + — 60 

4t 16 v ' 



11 Throughout this section we are referring to the diagonal geometry only, i.e., the transfer direction is diagonal 
with respect to the medial lattice. 
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Figure 17: Effective central charge for the loop model on the antiferromagnetic critical line in the 
sector N odd (diagonal geometry). The solid line shows the expected result, c e g- = c — 24/i£> = 

(2-f)-24(^). 

showing that the contribution from the boson <p\ is as in the N even case, and thus it must be 
the non-compact boson fa that is twisted, i.e., sees antiperiodic boundary conditions. 

Notice that if the field fa is twisted in the TV" odd sector, the continuous part of the spectrum 
must disappear, as the twisted sector of a compact and non-compact bosons are identical. The 
spectrum for N odd thus ought to be considerably simpler than for N even, a fact well in 
agreement with numerical studies. 

We now turn to the loop / cluster Potts model. The effective central charge is found numerically 
(see Fig. 117(1 to be given by a field of dimension 

hn = ^ (61) 
To interpret the presence of the weight (|61() in general, we observe that it can be written 

1 t (\ 1\ 2 1 



= i6 + 4 U-ij -« (62) 

As argued earlier, the contribution comes from an antiperiodic sector for the non-compact 
boson 02- The remaining contribution can then be interpreted as the weight for an electric charge 
eo = \ — j in the theory. Recall that in the ./V even sector we had charges eo = \ , 1 — \ , ■ ■ ■ 
instead. We do not know the origin of this additional \ contribution to the charge. 12 

12 Since in the numerical study, the twist for the vertex model is exactly a = h, it means that for N odd, this 
twist translates into an electric charge eo for the <j>\ boson, and a Zi twist for the fa boson. Moreover, the present 
value of the charge would correspond, in the usual case, to giving to an oriented non-contractible loop, a weight 
±ie ±t7T ^ . This says something about the relation between the microscopic arrow degrees of freedom and the bosons 
(j>\ and 4>2- 
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Figure 18: Effective watermelon exponents, in the sector N odd, for the loop model along the 
antiferromagnetic critical line. Diagonal geometry. 

Effective watermelon exponents in the twisted vertex model (and in the loop model) are then 
given by measuring the gaps with respect to the odd-iV ground state. This yields 

s (L/2) 2 -l 4-t 
h L =h L - h Do = — + -j^r- (63) 

and is verified numerically in Fig. 1181 

For t = k + 2 an integer, the above results can be fitted in a parafermionic picture by 
introducing the disorder operators (with respect to charge conjugation symmetry), whose weights 
read [2HJ 

\k - 2 + (k - 2s) 2 ] k , „, 

h °-- • 0£s£ 2 (64) 

The largest of these weights is obtained for s = 0, and coincides with (|61|) . On the other hand, it 
is well-known that correlation functions in the disorder sectors have to be defined on a two-sheet 
Riemann sphere, implying in particular the presence of descendents on half-integer levels. We 
expect this feature to extend to t generic, and this is indeed well verified by the numerics for the 
first excitations in the odd-iV sector, both for the loop model and for the twisted vertex model. 
We spare the reader the corresponding figures. 

8 Other exponents of the antiferromagnetic Potts model 

We consider first the magnetic exponent. In the Berker-Kadanoff phase, it is obtained by putting 
an electric charge at infinity so as to cancel the weight of non-contractible loops (a = ^), and 

"16 4t • (65) 
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Figure 19: Magnetic exponent of the loop model along the non-physical selfdual line. Axial 
geometry. The solid line represents the result (|65|). 

The dimension in the untwisted boson theory is thus A^ K = The dimension on the antiferro- 
magnetic line, conjectured first in [2], is obtained by adding a contribution from the non-compact 
boson, so 

"« =l6-«— i6-=8-4T (66) 
This value is easy to understand following the discussion after formula (|48|) : the magnetic expo- 
nent corresponds indeed to a = | (a value which kills non-contractible loops, in relation with the 
term having T](E) = 1 in ©), for which the XXZ value is never valid as t > 2 always, while the 
correct value deduced from the parafermionic exponents gives c e fr = —1, in agreement with the 
conjectured value of hff. 

The results (|65[) - ([66|) for the magnetic exponent are verified numerically in Figs. I19H20I 
Next, we consider the thermal exponent. An increase of K moves the model into a massive 
phase, while a decrease moves it into the Berker-Kadanoff (BK) phase. Numerical evidence is 
that the BK and massive phases are not related analytically in any way, and that the largest 
eigenvalue in the masive phase becomes the smallest (or maybe one of the smallest, in the sense 
that it scales with the smallest) eigenvalue in the BK phase. As a result, the free energy per 
unit area for the Potts model at Q generic and the twisted vertex model exhibits a discontinuity 
of the first derivative at the antiferromagnetic critical point, a characteristic of first-order phase 
transitions. Nevertheless, we have amply argued that these models right at the critical point 
have algebraic decaying correlations. We are thus in a case of a first-order critical point with 
an exponent v = 1/2 all along the line for the twisted vertex model, and, generically, for the 
loop/cluster Potts model. 

On the other hand, these results do not apply to the loop/cluster Potts model in the case 
of Q a Beraha number, that is t = k + 2 integer. For such values of Q, the free energy does 
not exhibit a kink any longer, the eigenvalues crossing the ground state at the antiferromagnetic 
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Figure 20: Magnetic exponent of the loop model along the critical antiferromagnetic line. Axial 
geometry. The solid line represents the result (|66j) , 

critical point being cut-off from the partition function within the Berker-Kadanoff phase. The 
same is true for the RSOS version of the model. These different features are illustrated in Fig. 1211 
While the existence of the first-order critical point can be pretty much argued analytically 
for the twisted vertex model and the Potts model, it is natural to expect a similar transition 
in the untwisted vertex model, since for the vertex model, the free energy per unit area should 
behave "normally" and thus be independent of twists. This is confirmed by numerical study, as 
illustrated in Fig. 

At the Beraha numbers, and for the loop/cluster Potts model or the RSOS model, the thermal 
operator does not have dimension | any more. Rather, the dimension becomes h = = |^|, 
the dimension of the parafermion operator. It is easiest to find out what happens for the RSOS 
version of the model, which has the same largest eigenvalue as the loop/cluster model. In this case, 
recall that the antiferromagnetic critical point is in the universality class of the Zk theory, and the 
flow in the Berker-Kadanoff phase coincides with the flow predicted by Fateev and Zamolodchikov 
|29j who studied perturbed parafermionic theories with the following hamiltonian: 

H = H k + \ (e ie / k iPi{j + e-^V^) (67) 

(A positive by convention). These authors argued that for 9 = tt the theory flows, for k odd, to 
the minimal model with central charge c = 1 — n^jn^m while it remains massive for k even. 

Numerical evidence for these flows is presented in Figs. l23H25| where we also give phenomeno- 
logical interpretations of other relevant features in the RSOS model phase diagrams. 

Note meanwhile that if t is rational, we can still define an RSOS model, but its ground state 
energy should be the same as the one of the vertex model — that is, exhibit the first-order critical 
point, in sharp contrast with what happens for t integer. This is illustrated in Fig. [23 for t = | . 

Finally the question remains of how to describe the first-order critical point within the 
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Figure 21: The dashed (resp. solid) curves represent the free energy of the RSOS model (resp. 
the loop model) at t = 5 for various widths N. The two free energies coincide perfectly (and 
hence, also with that of the twisted vertex model) for x = Q^ 1 ^ 2 (e K — 1) large enough (whence 
the dashed curves are 'hidden' by the solid ones). However, for smaller x, the free energy of the 
loop model experiences a sharp singularity with a discontinuity of the derivative. As N increases, 
this singularity moves towards the bulk antiferromagnetic transition temperature, here shown as 
a vertical line. Within the BK phase, the free energies of the RSOS and the loop models are 
different, with a jump whose magnitude increases as one goes deeper within this phase. Diagonal 
geometry. 
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Figure 22: The dotted (resp. full) curves represent the free energy of the untwisted (resp. the 
twisted) vertex model at t = 5 for various widths N. The dashed curves are for the RSOS 
model, same as in Fig. 1211 Only the twisted vertex model exhibits singular behavior for finite 
N. However, the free energies of the twisted and untwisted vertex model behave similarly in 
the thermodynamic limit, with a first-order critical point on the antiferromagnetic critical line. 
Diagonal geometry. 
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Figure 23: Parafermion to minimal model flow in the RSOS model with t = 5. The parameter 
x = Q^ 1 ^ 2 (e K — 1). The vertical lines show the non-physical selfdual line (x = — 1) and the 
antiferromagnetic critical point (x ~ —0.51). The horizontal lines are at c = 4/5 and c = 7/10, 
the central charges of the relevant parafermion and minimal models. The data indicates the 
existence of a new repulsive fixed point at x ~ —0.83 in the parafermion model universality class. 
The minimal model fixed point is situated somewhere between this and the antiferromagnetic 
point, and its attractive nature leads to a plateau in c in the range —0.83 < x < —0.51. The data 
further reveal a c = fixed point, presumably repulsive, at x ~ —0.93. If true, consistency of the 
RG flow diagram would necessitate another attractive fixed point; the data is indeed consistent 
with such a point at x ~ —0.87, again in the minimal model universality class. The interpretation 
of the data in the range —1 < x < —0.93 is somewhat uncertain. Right at x = —1, the data are 
again well-behaved and converge to c = —0.385 ± 0.001. We have no explanation for this value 
at present. Diagonal geometry. 
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Figure 24: RSOS model phase diagram for t = 6. The physics is different — and simpler — than 
in the case t = 5 shown in Fig. 1231 The c = 1 parafermion fixed point at the antiferromagnetic 
transition, x ~ —0.58, is still repulsive, but the nearest attractive fixed point is now the non- 
critical c = point at x ~ —0.87. Right at x = —1 the data converge slowly towards c = f .6±0.f , 
an unexplained value. The model is expected to be massive everywhere, except at the two vertical 
lines. Diagonal geometry. 
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Figure 25: RSOS model phase diagram for t = 7. The physics is similar to the case of t = 5 
shown in Fig. 1231 In particular, we observe again a parafermion (c = |) to minimal model (c = |) 
flow. The plateau at c = | develops more slowly in N, in the form of a 'shoulder'. Analyzing the 
behavior at the left edge of that plateau would call for larger system sizes than those presented 
here. Diagonal geometry. 
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Figure 26: The free energy of the RSOS model for t = | also develops a first-order singularity in 
the thermodynamic limit. Diagonal geometry. 
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Coulomb gas. We have not found the complete answer to this question, but a probable sce- 
nario is that we are dealing with the perturbation described by Fateev f2Q] > 



H = H + Acos(/3i$i)cosh(/3 2 $ 2 ) (68) 

where Hq is the hamiltonian for a pair of free bosons $1 and $2 as m section 6, and 

(3j - I3 2 2 = 4vr (69) 

Indeed, consider first the case where the boson $1 is twisted so that the central charge is c = 2— |, 

and set g = §■ By standard Coulomb gas calculations one finds then the dimension of the 

perturbation to be h = ^ 2 ± | = 0, 1, i.e., a mix of a perturbation with dimension zero and 
one with dimension unity. Meanwhile, the theory is integrable, and admits non-local conserved 
currents, one of them being 

(4i7T \ 
±-^ 2 J (70) 

with dimension h = 1 — = I — = frf > the dimension of the parafermionic current. [The 
other current meanwhile is 

4 = ^ 2 exp(±^x) (71) 

with dimensions h = 1 and h = ^4^-] A priori, the integrable perturbation involves a mixture of 
all these operators. What must happen is that for t integer, the quantum group truncation leaves 
only the contribution from J±, leading to the parafermionic perturbation. Otherwise, the two 
relevant operators have dimension zero and |5§ , presumably in agreement with our observation 
of a first-order critical point for the generic Potts model. 

Now in the untwisted theory, the perturbation in (|68|) has dimension h = ^, the one of J± is as 
before, and the one of J'± is h = One might be concerned by the absence of a dimension zero 
operator then, since we have argued that the untwisted vertex model had a first-order critical 
point as well. The point is that the coupling constant A in (|H5|) must be identified with the 
lattice coupling K — K c within the twisted theory, as multiplying the field of dimension zero in 
the twisted theory. Let us call it A tw for clarity. Then the physical dimension of A tw is [£] -2 . 
In the untwisted theory, the perturbation has physical dimension [L] -2 ' 1 = [X] -1 , so the coupling 
there must also have dimension i.e., be proportional to the square root of A tw , or the square 

root of K — K c . A singularity of the free energy such as \K — K c \ turns into | A tw | i.e., A 2 ntw from 
dimensional counting, corresponding, within the untwisted theory, presumably to A untw |A untw |, 
i.e., a thermal exponent a = 0, and thus an exponent v = ^, a perturbation of dimension h = \ 
indeed. 

Making this more precise is left for future work: it seems clear in any case that the variable 
K is not very natural to study the perturbations of the untwisted theory. 

A last observation we can make is that, since the parafermionic theories for t integer flow 
to minimal theories in the BK phase, and since the latter are derived from a Coulomb gas with 
coupling g = ^p, it is most likely that excitations for Q generic around the "false ground state"— 
i.e., the state that crosses from the higk K phase, and that would still be the ground state of 
the Potts model for Q a Beraha number — are described by this Coulomb gas, at least in the 
vicinity of the antiferromagnetic transition. Excitations around the true ground state meanwhile 
are described by a Coulomb gas with g = %, so both manifolds are in fact present within the BK 
phase! 
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9 Conclusion 



The least one can say is that the properties of the antiferromagnetic Potts model are extremely 
complex, even for the simple case of the square lattice. There are definitely loose ends in our 
study, the most important one being probably our lack of understanding of the emergence of 
a non-compact bosonic degree of freedom. One might understand it at a very qualitative level 
by considering the block interactions as in Fig. ^ Summing the arrows on either of the pairs 
of incoming or outgoing legs defines dual height variables which are described by the boson fa. 
The states for which the two arrows have opposite directions do not change the height fa, and 
draw loops on the lattice, which can intersect at vertices. These loops bear some resemblance 
to the loops of the Goldstone phases in 0(n) models, for which it is known that the continuum 
limit is essentially a collection of non-compact bosons and symplectic fermions. It is thus not 
so far fetched to expect that the loops in the antiferromagnetic Potts model are described by a 
non-compact boson, which would mean roughly that they are Brownian at large scales. It would 
of course be most interesting to elaborate this picture further. 

A more straightforward direction of study would be to consider, when t is rational, the RSOS 
versions on the antiferromagnetic critical line. Since when t = k + 2 is an integer, they are in the 
universality class of parafermions, identical with the SU(2)k/U(l) coset model, it is tempting 
to speculate that they might be related to SU(2) with a rational level. This will be discussed 
elsewhere. 

The appearance of a line of first-order critical points is intriguing. Such points are not 
entirely unheard of. In early works, they were exhibited for instance in the one-dimensional q- 
state classical clock model with a topological term |31| . Another case which is related is the point 
H = in the high-temperature loop version of the low-tempeature Ising model (i.e., the 0{n = 1) 
dense loop model [22] )• R- ecai l indeed that the low-temperature Ising model is a massive theory 
as far as the local spin degrees of freedom are concerned, but nevertheless presents interesting 
conformal properties when one reformulates it as a high-temperature loop expansion. These 
properties are the same as the ones of the dense 0(n = 1) loop model, and are described by a 
conformal field theory with vanishing central charge, and, more importantly, a field of vanishing 
dimension — the spin operator. Turning on a magnetic field leads to a first-order singularity of the 
free energy, that is, the model is also at a first-order transition point where the two manifolds of 
possible spontaneous magnetization m = ±m sp intersect. On the other hand, it is not clear how 
to see the difference between positive and negative magnetic fields within the loop formulation, 
since the field is conjugate to the number of loop end points, and that number is necessarily even. 
Therefore, it is not clear how to actually observe for instance the first-order singularity of the 
free energy. 

Interestingly, first-order critical points have also been independently proposed in a series of 
papers by Pruisken and collaborators revisiting the quantum Hall effect transition. These authors 
have argued for instance j3Hl that in the large- N CP N model, there are bulk massless degrees of 
freedom at the first-order phase transition point. If true, these massless degrees of freedom could 
maybe be understood in the Q-state Potts model incarnation as geometrical degrees of freedom, 
probably non-local in terms of the original Potts spins. 

Finally, the reader might wonder whether the antiferromagnetic transition discussed here 
is particular to the square lattice. We believe on the contrary that it is common to any two- 
dimensional lattice. Indeed, it is well established that the properties along the ferromagnetic 
transition line are universal, for example by invoking the standard one-boson Coulomb gas con- 
struction. The existence of another line of transitions for which the thermal operator is irrelevant 
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Figure 27: Central charge for the cluster model on the triangular lattice. The transfer direction 
is perpendicular to a line of N Potts spins. All fits are for N which have the same remainder 
modulo 3 (otherwise deviations from the expected result c = 2 — | will appear close to Q = 
and Q = 4). 

(which we have called, for the square lattice, the non-physical selfdual line) follows by analytic 
continuation in the Coulomb gas coupling constant, and leads immediately to the existence of 
a Berker-Kadanoff phase, in the region where e K < 1. To isolate this phase from the trivial 
fixed points at zero and infinite temperature, we need (at least) on either side a further line of 
repulsive fixed points. These lines must necessarily be the loci of level crossings, and it is thus 
not far fetched to expect that the rest of the physics of the first-order critical points (at generic 
Q) will follow. Moreover, the singularities at the Beraha numbers follow from the fact that the 
vertex model transfer matrix can be built in terms of Temperley-Lieb generators and that this 
commutes with the generators of the quantum group, two very universal ingredients. 

To test this heuristic argument, let us consider the case of the triangular lattice. Setting 
v = e K — 1, the Potts model is integrable when + 3v 2 = Q [HI]- When Q E [0,4] this 
describes three branches. The upper one is the ferromagnetic transition, of course, and by analytic 
continuation the middle one should then govern the Berker-Kadanoff phase. The remaining lower 
branch is thus a natural candidate for one of the two antiferromagnetic transition lines. Numerics 
for the central charge of the cluster model along the lower branch is shown in Fig.|27]and is indeed 
well described by c = 2 — |. The same results of course hold by duality for the model on the 
hexagonal lattice. 

Acknowledgments: J. L. Jacobsen thanks J.-F. Richard, J. Salas, and A. Sokal for related collab- 
orations, and SPhT for their warm hospitality. 
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A Block spins and the Q = case 



We use for instance the osp(2/2) R matrix given in [SSI- ^ n the latter reference, a more general 
expression is given for the g-deformed case. We take the limit g — > 1. At the same time we scale 



the spectral parameter x = e t/yv , q 



R = R 



e* 7 , with 7 
2 + v 



Po-V2 



0. We then obtain the expression 
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The R matrix acts in the tensor product of two fundamental, four-dimensional representations. 
The basis states in the latter are the bosonic states |1), |4) and the fermionic states |2), |3). The 
projectors follows: 

(i) One-dimensional blocks: 

Pi = 1, ^0 
Pi = 0, P 

(ii) Two-dimensional blocks: 
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(|1)®|4),|2)®|3),|3)®|2),|4)®|1)). 

Consider now a different problem where spins 1/2 interact with the well-known six- vertex 
model R matrix with anisotropy parameter 7. For the value u of the spectral parameter the 
Boltzmann weights are represented in Fig. I28| with 
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w 5 



W2 
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A sin 7 _ iu 
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s m 7 e - iu 



(78) 
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where A is a parameter at our disposal (a gauge parameter), and for A = 1 the R matrix commutes 
with su(2)q, with q = e* 7 . 



W, 



Figure 28: Weights in the six-vertex model. 



We now consider the following situation. Take a square lattice, and with every horizontal 
line associate the parameters h, h + ^ in alternance. With the vertical lines, associate similarly 
the parameters v,v + ^ in alternance (see Fig. I29j) . The weights at the vertices are given by 
the foregoing six-vertex formula with u being the difference of the horizontal parameter and the 
vertical parameter for the two intersecting lines. Since the Boltzmann weights are invariant under 
a shift of u by ir, we end up with two kinds of vertices, those for which u = h — v and those for 
which u = h — v ± ^ . 

We now consider interaction blocks where a pair of horizontal lines meets a pair of vertical 
lines. We use these blocks to define a new interaction matrix, where each link can carry four 
states, made out of the two spins 1/2 of the block. In the limit q — > i (7 — > ^), the interaction 
becomes trivial at fixed u, but if we let at the same time u approach zero as 



7 
u 



TT 

2 

euu 



+ e 



(79) 



we claim that we obtain, after inserting appropriate gauge factors, the osp(2/2) R matrix, with 
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Figure 29: Staggered vertex model, with alternating parameters on horizontal and vertical lines, 
the following correspondences 

= |1> 

I + -) + *!-+> = |2) 

| + ->-i|-+> |3) 

I--) = |4> (80) 

together with v = —2w. 

The check is laborious. We give here a simple example based on the fact that 



R osp (\l) |4» = |1> ® |4> + |4> |1) 



+^^2 (|2> |3> + |3) si |2» (81) 

This leads among others to the two vertices represented in Fig. I3UI These vertices are obtained 
in turn as a sum over configurations of the spins 1/2. For the first, one has the weight 

sin u \ 2 / cos u \ 2 w 2 



sin(7 + u) ) \cos(7 + n)/ (w + l) 2 



and for the second 



sin 2 7 2iu sin 2 7 sin 2 u 2 iu -, ^ 2 1 + 2w 

-e n- ; — -e 



sin (7 + it) sin (7 + n) cos 2 (7 + n) (to + l) 2 (1 + w) 2 

in agreement with the matrix elements from R sp- 

It is interesting to notice that as q — > i, the tensor product of two fundamental representations 
of su(2) q , which decomposes as the sum of a one- and a three-dimensional representation like for 
q = 1 when q is generic, becomes indecomposable, and has the structure shown in Fig. I3H where 
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(v-iy 



= -4 



(v-1) 



(v-2Y 



Figure 30: Two examples of the equivalence between R matrix elements. Diagrams on the left are 
osp vertices together with their Boltzmann weight, and they are obtained from the corresponding 
"microscopic" configurations on the right. 
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Figure 31: Indecomposable representation 1/2® 2 for q = i. 



\ui) = e J7r / 4 | H — } + e i7r / 4 | — h) and [^2) = e l7f//4 | H — ) + e J?r / 4 | — h), and the arrows indicate 
the action of the su(2) q generators: 

S~ = q° z O a~ + a' <g> 

5+ = ® (j + + cj+ ® (82) 
It turns out that, if one introduces the other generators 

5~ = q~ aZ a~ + a~ ® 

S + = q- a * ® a + + a + ® g aZ (83) 

&t q = i, together with the limits lim g ^j g+g-i » one obtains the four-dimensional representation 
ofoap(2/2). 
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